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ABSTRACT 

We examine two extreme models for the build-up of the stellar component of luminous elliptical 
galaxies. In one case, we assume the build-up of stars is dissipational, with centrally accreted gas 
radiating away its orbital and thermal energy; the dark matter halo will undergo adiabatic contraction 
and the central dark matter density profile will steepen. For the second model, we assume the central 
galaxy is assembled by a series of dissipationless mergers of stellar clumps that have formed far from 
the nascent galaxy. In order to be accreted, these clumps lose their orbital energy to the dark matter 
halo via dynamical friction, thereby heating the central dark matter and smoothing the dark matter 
density cusp. The central dark matter density profiles differ drastically between these models. For 
the isolated elliptical galaxy, NGC 4494, the central dark matter densities follow the power-laws r~°' 2 
and r -17 for the dissipational and dissipationless models, respectively. By matching the dissipational 
and dissipationless models to observations of the stellar component of elliptical galaxies, we examine 
the relative contributions of dissipational and dissipationless mergers to the formation of elliptical 
galaxies and look for observational tests that will distinguish between these models. Comparisons to 
strong lensing brightest cluster galaxies yield median (M*/£)b ratios of 2.1 ± 0.8 and 5.2 ± 1.7 at 
z k 0.39 for the dissipational and dissipationless models, respectively. For NGC 4494, the best-fit 
dissipational and dissipationless models have (M*/L)b = 2.97 and 3.96. Comparisons to expected 
stellar mass-to-light ratios from passive evolution and population syntheses appear to rule out a purely 
dissipational formation mechanism for the central stellar regions of giant elliptical galaxies. 
Subject headings: galaxies: formation — galaxies: elliptical and lenticular, cD — dark matter 



> 

in 
oo 
in 
o 

(N 
O 

o 



X 



1. INTRODUCTION 

How is the stellar content of galaxies assembled? Two 
extreme, contrasting models for the assembly of stars 
in galaxies have been considered over the years with no 
conclusive evidence as yet as to which mode dominates 
in which systems at which times. At one extreme, we 
can treat the process as totally dissipational with regard 
to energy; gas flows in from the virial radius, radiating 
away the kinetic and thermal energy it acquires while de- 
scending into the deep potential well of the dark matter 
halo. Once in place, the gas is transformed into stars 'in 
situ'-in approximately the regions in which we see stars 
in fully-formed galaxies today. The other extreme model 
postulates that stars are formed in smaller stellar systems 
far outside the effective radius of the ultimate galaxy. 
From there, they lose orbital energy via dynamical fric- 
tion and lose their potential energy only by heating or 
expelling other matter. These stars could be considered 
'accreted,' whether by minor or major mergers. The bal- 
ance between these two processes of galaxy formation 
(both of which surely occur) is unknown; determining 
that balance should help us to unravel some apparent 
paradoxes in the standard ACDM model of cosmology. 

Although the ACDM model of cosmology has been 
very successful explaining large scale observations, such 
as the cosmic microwave background and the large scale 
structure of galaxies, it has not enjoyed as much suc- 
cess on smaller, galactic scales. One notable discrep- 
ancy between simply-modeled ACDM predictions and 
observations is the difference between the predicted and 
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observed dark matter content in the centers of galax- 
ies. Many N-body simulations of dark matter parti- 
cles have been performed, and they show that cold dark 
matter will co llapse into self-similar halos with a central 
densit y cusp (iFukushige fe Makinolll997t iNavarro et al.l 
1997t iMoore et al.l 119991: iDiemand et all 120051: but see 



Springel et alj2008[) . The steepness of the cusp may vary 



from r 1 in the case of the NFW profile (Navarro et al. 
11997ft to as steep as r" 15 (IMoore et all 11999ft. Other 
studies (|Subramanian et al.l 120001 : iRicottil 12003: have 
shown a range of central slope indices, 1 < a < 1.5 for 
r~ a , depending on time, mass, and environment. 

Cold dark matter simulations uniformly predict that 
dark matter halos should be universally cuspy, but ob- 
servations have yet unambiguously to find these cusps. 
Indeed, observations of gravitational lensing, stellar ve- 
locity dispersions, and gas dynamics suggest inner dark 
matter density profiles are cored (r a , a > — 1) instead of 
cuspy dFlores fe Primacklll994HRomanowskv et al.ll200l 
[ Swaters et al Jl2003t iGentile et al.ll2004HSand et al.ll2004l: 
i Simon et al.1 120051: IC appell ari et aD 12006 1 : iGilmore et al. | 
20071: I Gentile et all 120071: Ide Blok et al.lbooa lOh et al l 
20081: iNapolitano et al.l l2009t but see iRhee et al.l 12004 



Spekkens et al.l 120051: iValenzuela et all 12007ft . For the 
Milky Way, it has been shown that within the errors of 
the microlensing observations, stars alone can more than 
account for the total mass density of the galaxy, leaving 
little room for a dark matter component in the center 
(jBinnev fe Evansll200lft . 

Acceptance of both the dark matter simulations and 
the observations of dark matter density profiles leaves 
us with an apparent discrepancy. One way to solve this 
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discrepancy is to abandon ACDM. For example, the self- 
interactions of warm dark matter will sm o oth central 
densi ty cusps ()Spergel fc Steinhardtl 1200 0; Bod e et all 
l2001f l . Another solution to this discrepancy may lie in a 
misunderstanding of galaxy formation and the assembly 
of stellar material at the center of dark matter halos. The 
addition of baryons to dark matter halos can have a pro- 
found effect on the dark matter profile of a galaxy. One 
such effect is the adiabatic contraction of dark matter by 
the slow addition of baryons to the center of the poten- 
tial well (|Blumenthal et a flll986). Adiabatic contraction 
has been obse rved in simulations with both dark matter 
and baryons (|Nayarro fc Bend [19911: IN avarro fc White! 
1994; Jcs seit et al. 2002tlAbadi et al. 120031: IGnedin et al.l 
20041) . This process conserves the adiabatic invariants of 
the dark matter orbits. With regard to energy, how- 
ever, it is a dissipational process, as the orbital energy 
of the infalling baryonic material is radiated away and 
lost from the system. Because the dark matter den- 
sity increases in the center of the galaxy under adia- 
batic contraction, the discrepancy between simulations 
and observations is only made worse, in the sense that 
the inner slope would be steeper than the NFW slope. 
For example, if the ultimate mass profile is "isothermal," 
{ptotir) °c r~ 2 ), then a dark matter profile with an ini- 
tial slope index 1 < a < 1.5, {r~ a ), would have a post- 
adiabatic-contraction index of 1.67 < a < 1.8. 

Processes which lower the dark matter density in the 
central regions of galaxies have also been explored. These 
processes are dissipationless; energy from stellar baryons 
is transferred to the dark matter, heating it and low- 
ering the central dark matter density. Such processes 
include interactions of the dark matter with a stellar bar 
([Weinberg fc Katdl2002t iHollev-Bockelmann et~aTll2005l : 
McMi llan fc Dehnenll2005fl , baryon energy feedback from 
AGN ( iPeirani et al.l 120081 ) , decay of binary black hole or - 
bits after galaxies merge ([Milosavlic vic fc Merrittj f2001), 
scattering of dark matter part icles by gravitational 
heati ng from infalling subhalos ([Ma fc Bovlan-Kolchinl 
200 1 ;. and dynamical friction of stellar/dark mat- 
ter clumps against the s mooth back ground dark 
matter halo (lEl-Zant et al.1 I2001L 120041: iTonini et al 
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20091) . In order to be successful, these methods must re- 
tain the strong central stellar concentration observed in 
la rge galaxies. 

lEl-Zant et al.l (|2001l ) propose erasing the dark mat- 
ter cusp via dynamical friction of incoming stellar/dark 
matter clumps aga i nst th e dark matter background. 
iRomano-Diaz et al.l ([2008D tested this hypothesis with 
dark matter/baryon N-body/SPH simulations and found 
that the introduction of baryons c an flatten the dark 
matte r cores in the inner 3 kpc. IRomano-Diaz et al.1 
(2008) claim that this cusp-flattening is not seen in other 
baryon/DM simulations because of a low numerical res- 
olution and a focus on ea rly, dissipational galaxy forma- 
tion. Recent simulations (iNaab et al.ll2007t lAbadi et al.l 
[2001 Uohansson et~aT1 2009) including both dark mat- 
ter and baryons have shown similar departures from 
the standard adiabatic contraction model. In their cos- 
mological simulations f or building up elliptical galaxies, 
Johansson et al. (2009) note that their results show rea- 
sonably low dark matter fractions in the inner 10 kpc 
and that the assembly of their elliptical galaxies at late 



times is dominated by accretion of stellar lumps, not 
gas, which as noted above, tend to reduce the dark mat- 
ter concentration. Similarly, accretion at late times has 
been invoked to erase the dark matter cusp formed by 
early adiabatic contraction, bring ing dark matter halos 
to a universal (e .g. NFW) profile ([Loeb fc Peeblesll2003l : 
IGao et al.ll200l . 

Dynamical friction of stellar lumps against the dark 
matter halo represents a dissipationless method of stellar 
build-up in which the orbital energy of the infalling stars 
is transferred to the dark matter, thereby heating it and 
driving it out of the central parts of the galaxies. This 
is in direct contrast to the case of adiabatic contraction, 
which arises from a dissipational build-up of stellar ma- 
terial. Although both processes undoubtedly take place, 
the fully dissipational (adiabatic contraction) and fully 
dissipationless (dynamical friction) scenarios for galaxy 
formation represent the two extrema. Since these two 
processes have opposite effects on the central dark mat- 
ter content of galaxies, the balance between them will 
determine the present-day dark matter density profiles 
in galaxies. 

In this paper we explore the physics of galaxy assem- 
bly, presenting two toy models for the two extremes of 
galaxy formation constrained to produce the same ob- 
served final stellar distributions: one model for the fully 
dissipational build-up of stellar matter with star forma- 
tion occurring in-situ and one model for the fully dis- 
sipationless build-up of stellar matter with stars added 
by accretion. The models are described in $2] We fo- 
cus on the structure of giant elliptical galaxies whose 
light profiles are well-described by Sersic models. Taking 
the observed stellar profile of the galaxy as given, but 
leaving the stellar mass-to-light ratio, M*/L, as a free 
parameter, we model the assembly of the stellar compo- 
nent via the two extreme formation methods. Then, by 
comparing the properties of the galaxies formed by these 
two methods to observations, we can begin to determine 
which method of assembling stars dominates and help 
to resolve the discrepancies between the simulated and 
observed dark matter density profiles. 

Throughout this paper, we adopt the ACDM cosmol- 
ogy with H = 70 h 70 km S -1 , fibaryon/^matter = 0.17, 

^matter = 0.26, and ft A = 0.74. 

2. MODELS 

In this section, we describe the two toy models of 
galaxy formation used in this paper. The first assumes 
that the orbital energy of the in-falling baryons is de- 
posited entirely in the dark matter halo, while the second 
assumes energy is radiated away, leaving the galaxy-halo 
system. 

In both models, we assume that both the dark matter 
and the baryonic matter follow circular orbits. For sim- 
plicity, we also assume the initial conditions (before the 
formation of the galaxy) for both the dark matter and 
the baryons are NFW density profiles, p oc [(?"/rNFw)(l + 
f/^Fw) 2 ] -1 ! with the same concentrations. The ratio of 
baryon to dark matter density is equal to the universal 

fraction, which we aSSUme tO be fibaryon/^darkmatter = 

0.20. In reality, the baryons will have a broader distri- 
bution than the dark matter, because the baryons are 
coupled to the radiation background. However, the ad- 
ditional thermal energy of the baryons corresponds to a 
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velocity less than 10 km s~ , which is negligible for ha- 
los more massive than ~ 10 8 M . In this study, we re- 
strict ourselves to large halos, and can safely assume the 
baryons are initially distributed identically to the dark 
matter. We have verified this is a valid assumption by 
placing all the baryons initially outside the virial radius 
where they have negligible binding energy. This does not 
affect any of the results in this paper, altering the dark 
matter mass with 5 R c by a few parts in 10 3 . 

The final condition to which the models evolve is an 
elliptical galaxy with a predetermined stellar luminosity 
profile in the center of a dark matter halo. In both the 
dissipational and dissipationless models, we assume that 
the final galaxies contain no gas; all the baryons are ac- 
counted for in stars. This assumption is justified because 
the gas mass in elliptical galaxies today is usually only 
a few percent of the tota l baryon mass in these galaxies 
(|Georgakakis et al.ll200l . 

For both models, the dark matter and stellar matter 
profiles are discretized into spherical shells such that the 
radius of a shell is at least two orders of magnitude larger 
than its width, and each shell is taken to be of homoge- 
neous density. In order to form the galaxy, shells of stel- 
lar matter are moved inwards. Each time a shell of stellar 
matter passes through a shell of dark matter, the stars' 
orbital energy is either deposited in the dark matter shell 
(dissipationless model) causing the dark matter shell to 
move outwards, or the stars' orbital and thermal energy 
is radiated away causing the dark matter shell's average 
radius to decrease as it undergoes adiabatic contraction. 
Since truly spherical shells would not suffer from dynam- 
ical friction, we are in fact considering a shell composed 
of stellar clumps, all of which (at any given time) have 
the same energy and angular momentum per unit mass, 
but whose angular momentum vectors are oriented ran- 
domly, giving the shell a total angular momentum of zero. 
Additionally, we assume that the stellar clumps do not 
contain any dark matter, as they formed far from the 
center of the dark matter halo, where the density is low. 
In reality, these clumps will contain some dark matter, 
much of which will be stripped before the stellar clump 
merges with the galaxy. Any remaining dark matter will 
undergo dynamical friction against the background dark 
matter. Because the system is spherically symmetric, a 
given dark matter shell is only affected by matter that 
crosses through it, but not by the rearrangement of mat- 
ter interior or exterior to it. Therefore a procedure that 
describes a single shell crossing can be repeated many 
times for many shells, until the desired galaxy has been 
assembled. 

2.1. Dissipationless Model: Dynamical Friction 

In the dissipationless model, all the orbital energy lost 
by the stars in forming the galaxy is deposited in the 
dark matter halo. In this case, the galaxy is built up by 
'dry,' dissipationless mergers. 

We imagine constructing the Sersic profile of the stars 
in a shell by shell fashion, sequentially moving each 
shell of stars from a large radius to its final position, 
in such a way that two stellar shells never cross. There- 
fore, the problem can first be idealized as moving one 
infinitely thin stellar shell from a large radius i?i to a 
smaller, final radius Ran- As the stellar shell moves, it 
passes through a dark matter distribution with a mass 



distribution Afdark(f), to which the stellar shell gives 
its orbital energy. The dark matter can be subdivided 
into a series of mass shells of thickness ARk, such that 
J2k ARk = Rint — Rfin- Therefore, we only need to com- 
pute the effects of one stellar shell passing through one 
dark matter shell uniformly distributed between Rk and 
Rk+i, with thickness ARk, and then execute a double 
sum over the dark matter shells and all the stellar shells. 
A depiction of this process is shown in Figure [TJ In the 
figure, the sum over dark matter shells is indexed by k, 
while the sum over the stellar shells is indexed by i. 




Fig. 1. — A cartoon depiction of the process by which we move 
stellar shells in to form a galaxy. The thin stellar shells are shown 
in red and labeled by Si. The dark matter background is shown 
in blue and discretized into layers of uniform density and width 
Rk ~~ Rk+i- The stellar shell Si is being moved from to Rk+i 
as shown by the arrow. 



The total energy lost by a single stellar shell as it moves 
from Rk to Rk+i can be calculated by taking the differ- 
ence between the energy change of the stellar shell as 
it moves from a large initial radius, i?i n t, with zero po- 
tential and zero kinetic energy, to Rk+i and the energy 
change as it moves from i?j nt to Rk- The changes in ki- 
netic and potential energy per unit mass of the stellar 
shell as it moves from i?i nt to Rk are 



1 GM(Rk) 

= 2 ~R k ' and W 

. „, GM(R k ) f R ^ GdM dark (r) 



where M(r) includes both the dark matter and any stel- 
lar matter interior to r. By assumption, there is no stel- 
lar matter between i?j nt and Rk, so the integral in the 
potential energy depends only on Mdark- The change in 
total energy of the stellar shell as it moves from Rk to 



4 



Lackner & Ostriker 



R 



fc+i 



is thus 



These energies are given by 



AE k 



GM dark (i? fc ) GM dark (i? fe+1 ) 

2R k 2R k+1 
GM star (i? fc ) /_i 1_ 

GdM dark (r) 



(3) 



The stellar shell is made of up of small lumps of stars 
that will undergo dynamical friction against the dark 
matter background as they move; the energy lost by the 
stars as they move from R k to R k +\ must be deposited 
in the dark matter. We assume that the exchange of en- 
ergy between stars and dark matter is a local process. 
Therefore, all the energy of the stars is deposited in the 
dark matter initially orbiting between R k and R k +i- As 
the dark matter gains energy, it's orbit will expand, and 
the uniform density dark matter shell between R k and 
Rk+i will become wider. Since we are assuming the en- 
ergy exchange is local, the dark matter orbiting at the 
inner radius (R k +i) will not move. In reality, this dark 
matter is affected by both adiabatic contraction, since 
the stars add mass interior to Rk+i, and dynamical fric- 
tion, since the energy exchange is not purely local; this 
layer will move, but the movement will be second order in 
AR k = R k — Rk+i- Therefore, by conserving energy the 
only quantity which changes is thickness of the dark mat- 
ter shell the stars have moved through and consequently 
its mean radius. 

In order to calculate the new dark matter layer thick- 
ness, we repeat the procedure above, but this time also 
account for the energy of the dark matter shell both be- 
fore and after the stars move through it. As above, we 
assume a spherical dark matter shell of mass A/dm and 
uniform density, with an inner radius of R (correspond- 
ing to R k +\ above), and a thickness Ai? « R. Directly 
exterior to the dark matter shell is an infinitesimally thin 
stellar shell of mass AM*, at the radius R + AR. The 
kinetic+potential energy of the two-shell system can be 
broken into the self-interaction energy of the dark matter 
shell, Ed, the self- interaction energy of the stellar shell, 
E a , and the interaction energy of the two shells, E ds - To 
first order in Ai?, these are given by: 



E d 
E s 
Eds - 



GM DM ( _ 2 AR 
3 R 
AR 
~R 

1 - 



1 



4i? 
GAM, 2 
4i? 

GAM*M DM 
2R 



AR 



(4) 
(5) 
(6) 



The shells are embedded in a spherically symmetric 
galaxy which also contributes to the energy. We assume 
that the mass distribution of the galaxy remains fixed as 
the shells interact. This assumption is exactly true for 
the mass interior to the shells and the mass far outside 
the shells. We define M- mt to be the total mass (dark 
matter + baryons ) interior to R. This mass contributes 
kinetic+potential energy E- lnt to the shells while the mass 
external to R + AR contributes a potential energy E ext - 



E;„t — — - 



2R 



Mdm 
AM* [ 1 



1-**) 
2R J 

AR 



and 



E r 



-G (AM* + Mdm) 



dM(r) 



(7) 



(8) 



R+AR 



Therefore, the total initial energy of the two shells is 
given by the sum of equations (J4j)- ((HJ) - 

We now move the stellar shell through the dark matter 
shell until the stellar shell is orbiting at the radius R, 
as in the above example. This time, however, we will 
include the changes in the dark matter distribution in 
the energy difference. After the move, the dark matter 
shell will widen from Ai? to Ai?'. Because Ai?' is still 
small compared to R, the density of the dark matter layer 
is still uniform after the move. The final total energy of 
the stellar and dark matter shells after moving the stars 



is 



E-. 



G_ 
2R 



Km 



1 - 



2 Ai?' 



AM; 



AM*M DM ( 1 - 
M int (M DM I 1 



3 i? 
Ai?' 



(9) 



2i? 
Ai?' 
~2R 



AM, 



-G (AM* + Mdm) / 

J R 



dM(r) 



AE S 



R+AR. 



If we assume that the mass initially exterior to R + AR is 
unaffected by the dark matter expanding to i? + Ai?', the 
change in external energy, AE cxt , depends only on the 
mass initially between i?+ Ai? and i?+Ai?', and the dark 
matter that moves beyond i? + Ai?. This assumption is 
valid to first order in Ai?. Since both Ai? and Ai?' 
are small compared to R, we assume that mass between 
i? + Ai? and i? + Ai?' is of homogeneous density and a 
total mass M cxt . Therefore, 



AE^t — 



GM DM M ext Ai?' - Ai? 



4i? 



Ai?' 



(10) 



Substituting this into equation [5] and taking the differ- 
ence between the total energy of the two shells before 
and after the move yields 



AE = 



G 



(Ai?' - Ai?) M DM 2Ai 
3 ( M int + AM, 



DM 



(11) 



i? 

~AR' 



M„ 



3 ARAM., (M DM + 2M in t + AM*) 



This can be solved numerically for Ai?', keeping in mind 
that M cxt is a function of Ai?'. In the limit that M ex t — »• 
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0, equation [TT] can be solved analytically for AR': 

AR' _ l/2M int + 1/3A/dm + 3AM, 
AR l/2M in t + l/3M D M + 3/2AM» 1 ' 
(AM,/M DM )(M int + 1/2AM«) 
l/2M int + 1/3M DM + 3/2AM* ' 

In the limit that AAf* << Mdm the ratio of widths goes 
to 1 as expected. This procedure can be easily scaled up 
to a series of interleaved dark matter and stellar shells. 
As the stellar shells move inwards to form a galaxy, they 
expand each dark matter shell they cross, thereby slowly 
moving dark matter outward. 

In the example above, all the energy from the stars 
is deposited in the dark matter layer which the stars 
cross. In our numerical calculations, the stars deposit 
their energy in a set of layers surrounding the layer they 
cross. The width of this set is proportional to the cur- 
rent radius and the amount of energy deposited in each 
layer is proportional to the mass of that layer. This ap- 
proximates the wake created by infalling stellar mate- 
rial, which is responsible fo r the forces causing dynami- 
cal friction (| Weinberd ll986) . The size of the wake scales 
as GAM*/cr| m w (AM*/M int )i?. We assume that the 
stars are added by a series of minor mergers, in which 
each added stellar shell of mass AM* constitutes a minor 
merger. The mass ratio of such mergers is approximately 
(AM*/Mj n t) ~ 1/10, thus setting the size of the wake. 

In a more accurate treatment, there would also be a 
diffusive term; each dark matter shell would spread out in 
radius as it was on average moved outwards when passed 
by a stellar shell. 

In the calculations above, we assume no stellar shells 
cross each other. This is equivalent to assuming all the 
infalling material remains on spherical orbits. This is not 
true in reality, especially since galaxies are not spherical 
but often triaxial systems, which do not allow purely cir- 
cular orbits as assumed above. If we allowed for triaxial 
systems in our models and allowed for the accretion of 
material on radial orbits, the infalling stars would de- 
posit energy interior and exterior to their final mean or- 
bital radius. This would have an effect on the final dark 
matter density profile, but the effect would depend on 
the fraction of energy deposited interior and exterior to 
the final orbital radius. If all the energy is deposited in- 
terior to the final orbital radius, then more dark matter 
would be displaced from the center, leading to a lower 
central dark matter density. The opposite is true if the 
energy is deposited outside the final orbital radius. 

Additionally, if we relax the requirement that no stellar 
shells cross, we must take into account energy deposited 
in the stars, not just the dark matter. This will expand 
the stellar orbits, in the same way the dark matter or- 
bits are expanded, and lower the stellar density in the 
center of the galaxy in the same way the dark matter 
density is lowered. Therefore, in order to make a galaxy 
with a given stellar density, we would first have to make 
a more concentrated stellar system, and then add stellar 
clumps which would undergo dynamical friction against 
the highly concentrated stellar system, thereby lowering 
the central stellar density to the desired value. Indeed, 
observations have been made o f highly concentrated stel- 
lar systems at higher redshift ([van Dokkum et all [20081: 
ICappellari et al.ir2009D . In order to make highly concen- 



trated stellar systems, the galaxies would have to un- 
dergo an early period of dissipational formation. This 
would also increase the central dark matter density be- 
fore the onset of dissipationless formation, making the 
final dark matter density dependent on the relative im- 
portance of dissipational and dissipationless formation 
mechanisms. Simulations show that early type galaxy 
formation can be divided into two phases, an initial dis- 
sipational formation of a centrally concentrated system, 
followed by a ccretion via 'dry' mergers of additional ste l- 
lar material ([Naab et al.l 120071 120091 : iCook et all 120091 ) . 
If we allowed stellar shell crossing, we would have to take 
into account the two phase growth of galaxies and com- 
bine the dissipational and dissipationless models into a 
single model. This work attempts to determine the rel- 
ative importance of the dissipational and dissipationless 
formation mechanisms; we are only concerned with the 
two extreme formation models, which can be easily mod- 
eled assuming spherical galaxies and infalling material on 
circular orbits. The effects of triaxiality and radial or- 
bits would require combining the dissipational and dissi- 
pationless formation mechanisms, which is left for future 
work. 

Of course, dynamical friction on incoming stel- 
lar clumps is intrinsically a t hree-dimensional process 
(ITremaine & Weinberg] 119841: IPichon fc Auberti 120061: 
lAubert fc Pichonl 120071) . and the treatment in spherical 
shells above is not intended to accurately mimic the ac- 
tual assembly of a galaxy via dynamical friction. Rather, 
the adopted model is designed to be correct with respect 
to the total energy deposited in the dark matter. In the 
numerical calculations presented below, the total bind- 
ing energy is conserved. As the width of the stellar shells 
decreases, the calculations conserve energy to approxi- 
mately one part in 10 5 of the binding energy of the stars 
in the final galaxy. 

2.2. Dissipational Model: Adiabatic Contraction 

For the dissipational build-up of galaxies we present 
the following picture: baryons in the form of gas slowly 
fall into the centers of dark matter potential wells, where 
the gas condenses to form stars. In this case, the gas 
radiates away its orbital and thermal energy as it falls 
inwards, so the total energy of the system is not con- 
served. However, as long as the gas is accreted into the 
center of the galaxy on a timescale that is long compared 
to the local dynamical time, the adiabatic invariants of 
the dark matter orbits will be conserved. In our model 
for adiabatic cont r action , we follow the prescription of 
iBlumenthal et all (j!986|) and again assume circular or- 
bits for the dark matter and a spherically symmetric 
mass distribution. Instead of energy being conserved, 
the adiabatic invariants of the dark matter orbits are 
conserved. For periodic orbits, § pdq is an adiabatic in- 
variant, where q is a coordinate and p is its conjugate 
momentum. For a particle in a circular orbit at a radius 
r around a spherical mass distribution M(r), we take the 
conjugate momentum to be the angular momentum and 
its corresponding coordinate, the angular position. The 
adiabatic invariant is then: 

J 2 = (I y/M(rjrde) oc rM(r) . (13) 
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Therefore, if the mass interior to the orbit, M(r), in- 
creases, the orbital radius must decrease. Following the 
same set-up as in the previous model, we start with a 
dark matter shell of constant density and mass Mdm di- 
rectly interior to an infinitesimally thin shell of stars/gas 
of mass AM*. After the baryons move interior to the 
dark matter layer, the inner radius of the dark matter 
shell becomes R! = R M int /(M int + AM*) and the width 
of the layer becomes 

AR> = (M int + M DM ) (fl+Afl) _ , f . 

M int + M DM + AM* ' 1 ' 

Overall, the dark matter shell moves inwards and be- 
comes thicker or thinner depending on the mass ratios of 
the dark matter shell, the baryon shell, and the mass in- 
terior to both. As in the previous model, two shells only 
interact when their orbits cross, so this prescription for 
interchanging two shells can be scaled up to many shells. 

If the constraint of circular orbits is removed, the adi- 
abatic i nvariant is no l onger rM(r). Using N-body simu- 
lations, [GnedhTeFjjl] ((2003) show that M(f)r, where f is 
the orbit-averaged position, is a good proxy for the adia- 
batic invariant. If this quantity is conserved for isotropic 
orbits, the prescription for adiabatic contraction remains 
the same. However, the prescription above will overes- 
timate the amount of adiabatic contrac tion if the orbits 
are radially biased (|Gnedin et all 120041 ). 

2.3. Models Fit to Example of Massive Galaxy 

The difference in the dark matter density in the cen- 
tral regions of a galaxy for each model is clearly shown by 
the comparison of the velocity profiles of model galaxies. 
In the top panel of Figure [2j the circular velocity curves 
for both models are compared to the velocity curve pro- 
duced by the stars alone. The models shown in Figure 
[5] are taken to fit very massive galaxies and we have ad- 
justed M*/L such that the central velocity dispersions of 
both models are the same. The central velocity disper- 
sions are computed assuming the stars are on isotropic 
orbits. In both cases, the stars dominate in the very 
central regions, but the differences in dark matter signif- 
icantly affect both curves. The lower panel of Figure [5] 
shows the ratio of the dark matter to stellar matter pro- 
jected densities as a function of radius for both models. 
In reality, since galaxies certainly form by a combination 
of dissipational and dissipationless methods, the velocity 
curves and the dark matter to stellar matter density ra- 
tios will fall somewhere in between the two extreme toy 
models examined here. 

In order to compute the two models, we require a set 
of input parameters describing the initial matter distri- 
bution and the final stellar distribution. If we assume 
the initial distribution of both components fit a single 
NFW profile, then the two required input parameters are 
the total mass of the halo and the initial concentration. 
Assuming that the galaxy is spherically symmetric, the 
final stellar mass distribution is completely described by 
the surface brightness profile of the galaxy and a stellar 
mass-light ratio. The luminosity profile can be directly 
observed and M*/L can be determined from observations 
of the central velocity dispersions of galaxies. Since stars 
dominate the mass in the central regions, they also are 
the dominant contribution to the central velocity dis- 
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Fig. 2. — The circular velocity (top) and projected dark matter to 
stellar density ratio (bottom) for an extremely massive L = 6.0 L, 
(it-band) galaxy in a 3.84 X 10 14 Mq halo. The Sersic index is 
n = 6.80 and the half light radius is 45.79 kpc. The stellar mass- 
to light ratios of the models are set to the values shown in the 
lower panel of Figure \3\ such that the two models have the same 
central velocity dispersion (cro ~ 270 km s -1 ). The solid line is for 
the dissipational model. The dashed line is for the dissipationless 
model. In the upper plot, the upper (dissipationless) and lower 
(dissipational) dotted lines are the circular velocity curves due to 
the stars alone. 

persion. Thus, the observational input parameters for 
these models are the total luminosity of the galaxy, the 
luminosity profile (in the case of Sersic profiles, the nec- 
essary terms are the Sersic index, n, and the half-light 
(effective) radius, i? c ), and the central velocity disper- 
sion. The free parameters are the total halo mass, the 
initial NFW concentration, and the stellar mass-to-light 
ratio (which we assume to be constant throughout a given 
system). Therefore, for a given galaxy, the two models 
will by construction have the same luminosity and cen- 
tral velocity dispersion, but the total halo masses and the 
mass-to-light ratios will be different, and in some cases, 
outside the bounds set by other observational and model 
constraints. Table [I] gives the parameters for the dis- 
sipational and dissipationless models which best fit the 
strong lensing cluster MS-2I37-23 discussed in £ 13.41 

2.4. Minimum Radius for Galaxies Formed by 
Dissipationless Mergers 

In order for a galaxy to form by purely dissipationless 
processes, the incoming stars must deposit their orbital 
energy in the dark matter halo. Therefore, the dark mat- 
ter halo must initially have sufficient binding energy to 
give to the infalling stars. The stars in the final galaxy 
can have no more binding energy than the initial dark 
matter halo. This sets a lower limit on the size of a galaxy 
formed by purely dissipationless processes. A galaxy 
formed by dissipationless accretion would approach the 
minimum size and would have very little dark matter in 
the center, which is in agreement with observations men- 
tioned in SJTJ However, if dissipational processes play a 
dominant role in galaxy formation, there is no minimum 
size for galaxies, as the infalling baryons can dissipate all 
of their orbital and thermal energy. 
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Given a cluster mass, concentration, and final galaxy 
stellar mass, and Sersic index, a lower limit on the 
galaxy's effective radius, R c , can be obtained by set- 
ting the change in binding energy of the dark matter 
halo equal to change in binding energy of the stars and 
baryons in the cluster. For galaxies with L w 10 11 Lq 
in the if-band, the minimum effective radius for a dissi- 
pationlessly formed galaxy yields a reasonable minimum 
for the effective radiu s of observed early type galaxies. 
(|Bernardi et al.ll2003fl . 

We should note that there is another limiting radius 
for the formation of accreted halos. A satellite may be 
destroyed by tidal shocks before it reaches the energeti- 
cally allowed minimum radius, thereby depo siting stars 
in th e outer region of th e growing galaxy (iKormendvl 
[19771 IGnedin et allflflQfll : I Wetzel fc White! 120090 . This 
process is difficult to compute but gives a limiting radius 
comparable to the dynamical friction limits computed 
above. For the most massive systems (BCGs), the tidal 
shock limiting radius is more severe than the dynamical 
friction one, so we can expect that the dark matter to 
stellar matter ratio in such systems will be higher than 
in more moderate mass galaxies. 

3. COMPARISONS TO OBSERVATIONS OF BCGS 

One test of the dissipational and dissipationless models 
of galaxy formation is the build-up of brightest cluster 
galaxies (BCGs). BCGs offer a good comparison sam- 
ple for these spherically symmetric test models for sev- 
eral reasons. First, BCGs almost always sit at or nearly 
at the center of their host cluster. Therefore, they are 
also centered in a dark matter halo, so there are no con- 
tributions to the potential from an off-center halo, not 
included in our models. Furthermore, BCGs represent 
a uniform sample, so much so that they have been sug- 
gested as standard candles ([Postman fe Lauerlll995l ). al- 
lowing comparisons to be made to the entire population 
instead of individual galaxies. Finally, BCGs are thought 
to have been formed by a se ries of galaxy mergers dur- 
ing the build-up of clusters (lOstriker fe Hausmanl 1X9771 ; 
iNipoti et all 120041: iCoorav fc Milosavlievicll2005[ ). mak- 
ing them good candidate systems in which to observe 
dissipationless galaxy formation. 

3.1. Scaling Relations for Input Parameters 

In order to compare the models to observations, we 
normalize the models to if -band and r-band data of 
BCGs. Given the luminosity of a BCG and choosing 
a constant stellar mass-to-light ratio for the models, we 
use empirical rel ations to derive the c luster and the BCG 
properties. From lLm fc Mobirl (|2004| ). the cluster mass is 
related to the observed if-band BCG luminosity by 



L 



BCG 



10 n h 



(4.9 ±0.2) 



70 



A/, 



cluster 



10 14 Mr. 



0.26±0.04 



(15) 



The luminosities of BCGs are 6-10 times brighter than L* 
in the i f-band. For galaxie s typically in clusters, Mk* = 
-24.34 (jLin fc Mohrll2004D . including the BCG, or about 
1.16 x 10 11 Lq in the if-band. From the cluster mass, 
the cluster virial radius, r2oo, can be calculated assum- 
ing the critical density p cr = 1.36 hy x 10 11 M Q Mpc -3 . 
Simulations have shown that the dark matter halo con- 
centration, c = r-200/^NFW; scales approximately with 



the mass as (|Neto et al.ll2007ft 



4.67 



M, 



luster 



10 14 M, 



(16) 



Equations Qj)] and [TH] set the initial conditions for the 
models. The properties of the stellar component of 
the BCG can also be derived from Lbcg- We assume 
that the BCGs are well-modeled by a single Sersic pro- 
file (I(i?) ~ expffl 1 /"). ISCTsia[l968h . ig noring the ICL 
and o uter components of the BCG (see [Gonzale z et al.l 
2005). The two-dimensional surface brightness profile 
defined by Sersic can be deprojected numerically into 
a three-dimensional luminosity density profile, assuming 
the galaxy is spherically symmetric. This numerical de- 

rojection is well- appro ximated by the analytic formula 

Lima Neto et al.lll999ft : 

p,(r)oc(r/i? c ) 1 - 1 - 188 /( 2 ") +a22 /( 4 " 2 ) 

exp ((0.327 - 2n){r/R c ) 1 / n ) , (17) 



where n is the Sersic index (n = 4 for a de Vaucoulcurs 
profile), and R e is the half-light radius of the surface 
brightness profile. Observations show that the Sersic 
properties of BCGs are corr elated with the gal axy's 
luminosity. Using d ata from ILin fc Mohrl (|2004f) and 
iGraham et al.1 (|1996l ). the luminosity can be related to 
the half-light radius and the Sersic index by 

log R c = -10.30+ 1.01 log (jj^j and (18) 



n = 2.91ogi? c + l-98 : 



(19) 



whic h are in good agreement with scaling rela tions 
fromlVale fc Ostrikerl I (120081) . iBernardi et al.ll200l and 
iDesroches et al.l (|2007T ). 

Also modeled in each galaxy is a central supermas- 
sive black hole, which adds a minor correction to the 
velocity dispersion of the galaxy. The black hole mass 
is determined f rom the galaxy luminosity by the relation 
(|Graharnll2007T) : 



log 



Ah 



BH 



Mr. 



(20) 



(0.95 ±0.15) log 



J BCG 



10 1 



L 



(8.26 ±0.11) 



Thus, by supplying a galaxy luminosity and a stellar 
mass-to-light ratio, we can obtain all the other input pa- 
rameters needed to compare the dissipational and dissi- 
pationless models to observations of BCGs. 

3.2. The L-a Relation 

The innermost probe of the mass profile is the cen- 
tral velocity dispersion of a galax y. Elliptical galax- 
ies fa l l on the fundamenta l plane |jDiorgovski fc Davisl 
Il987t iDressler et al.l I1987T) and one projection of the 
plane is the Faber- Jackson relation: the relation be- 
tw een a galaxy's luminosi ty, L, and velocity dispersion, 
o~ (jFaber fc Jackson! 119711 . In the if-ba n d, th e Faber- 
Jackson relation observed by iPahre et al.1 (|1998l ) is 



M K = -10.35 ±0.55 log cr . 



(21) 
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For BCGs. lLauer et ai1(|2007t ) show that the velocity dis- 
persion saturates at high luminosities, leading to the re- 
lation 



M v = -2.5(6.5 ±1.3) log 



250 km s~ 



-22.45 ±0.18. 

(22) 

iDesroches et al.l (|2007[ ) find a similar relation. In order 
to compare to the L-a relation for BCGs, we normal- 
ize our models to observed BCGs using the scaling rela- 
tions described in ij3.ll The line-of-sight central velocity 
dispersions averaged over an aperture of 1.64 kpc are 
then calculated for both the dissipationless and dissipa- 
tional mode l s. The comparison to the L-a relation from 
lLauer et al.l (|2007| ) is shown in Figure El Although both 
the dissipational and dissipationless models are slightly 
steeper than the L-a relation found by Lauer, the dissi- 
pationless model has a slope that more closely matches 
the observed L-a for BCGs. To match each model to the 
observations, the stellar mass-to-light ratios can be ad- 
justed. For the dissipational models, the best fit M^/L 
in the if-band is 1.43, while for the dissipationless mod- 
els the best fit to the L — a relation is for M»/L = 2.40. 
These are equivalent to stellar mass-to- light ratios of 7.57 
and 12.71 in the U-band, assuming V — K = 3.31 for the 
BCG population. In t he if-band, the st ellar mass to 
light ratio measured by iCole et ail ()200lD is 0.73 for a 
Kennicut IMF and 1.32 for a Salpeter IMF. 

The stellar mass-to-light ratios derived for these mod- 
els are simply the mass in stars needed to reproduce the 
dynamics (in this case, the central velocity dispersion) di- 
vided by the total observed luminosity of the galaxy. Al- 
though the dissipationally formed galaxies were brighter 
at high redshift due to star formation, we are only con- 
cerned with the z k luminosity and dynamical state 
of the galaxy. This corresponds to the luminosity of 
the evolved population; therefore, we have implicitly in- 
cluded passive evolution in the dissipational model and 
do not need to passively evolve the M*/L values derived 
above. 

However, the mass-to-light ratios are sensitive to the 
empirical scaling relations. For example, if the rela- 
tion for -B e (LBCo) is replaced with that derived by 
IBernardi et al.l (|2007T ) for galaxies fit by de Vaucouleurs 
profiles (n = 4), the best fit stellar mass-to-light ratios 
become 0.98 and 1.76 for the dissipational and dissipa- 
tionless models respectively. 

Additionally, the above calculations rely on the scal- 
ing relation between the total halo mass and the BCG 
luminosity. Instead, we can assume that the cluster is 
built hierarchically out of galaxies formed at z w 2. 
At z = 2, the concentratio n of a dark matt er halo is 
a weak function of halo mass; iGao et ail (|2008l ) find that 



c oc Af ha ° 031 . The fraction of mass in a halo which will 
form stars is given by M^/M^io oc M^^ 6 (|Lin et al.l 
120031: iBode et al.l 12009). Using these relations for the 
concentration and stellar mass, the best-fit M*/L values 
become 0.69 and 1.55 for the dissipational and dissipa- 
tionless models respectively. These values are in better 
agreement with the measured values given above. 

It is not surprising that neither value for M*/L can be 
discarded based on the Faber- Jackson relation. In the 
central regions of the galaxies, both models are stellar- 
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Fig . 3. — The L-a relation in the X-band from Lauer et al. 
(2007) compared to the calculated relation for both the dissipa- 
tiona l and dissipa tionless models. The dotted line is the relation 
from lLauer et al.l l l2007f) . The top fi gure shows both models at 
{M f /L)x = 1.60. By adjusting the stellar mass-to-light ratio for 
each model separately, the model lines can be made to overlay the 
observed relation, as shown in the lower panel. 

dominated, as shown in the lower panel of Figure [2j 
Therefore, even though the dark matter density can dif- 
fer by more than a factor of 10, it only makes up <~ 10% 
of the total mass in the central regions, and thus does 
not determine the central dynamics. 

3.3. Microlensing Optical Depth 

Microlensing of quasars, whic h has been observed 
in multiply imag ed systems ([Wozniak et al.l 120001 : 
Wambsganss 2006, and references therein), in prin- 
ciple provides a probe of the mass function of mi- 
crolenses (MACHOS, stars, or dark matter substruc- 
ture), and the density of the se microlenses relative to 
a smooth background density (ISchechter fc W ambsganss 
[20021: iDobTer et al.H2007t IPoolev et al.ll2009l) . Searches in 
the Sloan Digit al Sky Survey have found ~ 220 strongly 
lensed quasars (|Inada et al.l 120081 ). which are lensed by 
individual galaxies or entire clusters. In the following, we 
use the best-fit models for BCGs as an example to show 
the expected differences in microlensing results between 
the dissipational and dissipationless formation models. 

The difference in stellar mass-to-light ratios between 
the dissipational and dissipationless models leads to dif- 
ferences in the microlensing optical depth. The optical 
depth, r, is proportional to the number density of lenses, 
stars in this case, times the Einstein radius, #e, of each 
lens. Assuming that the distance between lens and source 
is large compared to the size of the ga laxy, the microlens- 
ing optical depth is ()Paczv nski 198g): 



r = XL 



4ttG AsA 



A 



(23) 



where is the projected stellar density and A are the 
angular diameter distances to the lens, to the source, 
and between the lens and the source. If (A S A) / A — D 
and the Sersic index of the lens galaxy is assumed to be 
4.0, then the microlensing optical depth at the half-light 
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radius is 



r=1.33 x IQ- 2 (M*/L) 



L 



BCG 



Hr- 



IO kpc 



K V 10 11 L 

2 , D 



(24) 



0.5 Gpc 



For the dissipational and dissipationless models that best 
fit BCGs, the ratio of the microlensing optical depths 
equals the ratio of M*/Z, or 1.43/2.40 = 0.60. In most 
cases, the microlensing optical depth at the position of 
the image is of order unity. Occasionally, individual mi- 
crolensing events can be observed. 

Additionally, the relative density of smoothly dis- 
tributed m atter (dark matter) to microlenses (stars) can 
be probed (ISchechter fe Wambsg anss 20021: iDobler et all 
120071 : iPoolev et all 12009( 1. Using the dissipational and 
dissipationless models shown in Figure [2] as an example 
(L = 6.0 L„), the fraction of dark matter to total mat- 
ter along a line-of-sight is 0.81 at 0.1 R c and 0.98 at 1.0 
R e for the dissipational model. For the dissipationless 
model, the ratios at 0.1 R e and 1.0 R e are 0.48 and 0.97, 
respectively. These large differences should be measur- 
able in microlensing studies of multiply-imaged quasars. 

3.4. Strong Lensing 

Strong lensing measurements provide a clean method 
of probing the total mass distrib ution of BCGs and their 
host clusters. [Sand et al. (2004) present observations of 
radial and tangential arcs for six clusters acting as lenses. 
The positions of the radial and tangential arcs are given 
by the solutions to 



= 1 



= 1 



d M pro j(i? rad ) 
dR 7ri? rad 

-^proj(-Rtan) 



and 



(25) 



where M pro j (R) is the projected mass interior to R scaled 
by the critical surface density, 



AttG Di s Di 



(26) 



Together, the radial and tangential lenses constrain 
the slope of the density profile and its normalization. 
ISand et al.l §004) use the lensing information as well as 
the velocity dispersion profile of the BCG to create den- 
sity models for the stars and the dark matter in each 
lens. They find that the mean inner dark matter density 
profile for six lensing clusters is r -°- 52 ±°- 3 j significantly 
shallower than the NFW profile. 

We repeat the analysis of ISand et al.1 (|2004| ). fitting 
our dissipational and dissipationless models for the dark 
matter profile to the three clusters with both radial and 
tangential arcs. As in the previous section, we assume 
that each BCG in the center of the cluster is built hierar- 
chically, either a series of purely dissipationless mergers 
of smaller stellar systems or by the dissipational accretion 
of gaseous streams which lead to in situ star formation. 
The dissipationless model for formation will yield a lower 
dark matter density in the center of the cluster, while the 
dissipational model for formation will lead to adiabatic 



contraction of the dark matter. At z ~ 2, both dissipa- 
tional and dissipationless formation mechanisms will be 
important, but the ratio between the two mechanisms 
is unknown, and by comparing data to the purely dissi- 
pationally and dissipationlessly formed BCGs, we hope 
to constrain how much each mechanism contributes to 
galaxy formation. 

For the dissipational and dissipationless models, the 
fixed input parameters for the models are the BCG lu- 
minosity, half-light radius, and Sersic index (n = 4). The 
free parameters are the total cluster mass, M c i uster , the 
dark halo concentration, c, and M*/L. By randomly se- 
lecting these input parameters from a reasonable range, 
we can find both dissipational and dissipationless mod- 
els that are within 1 and 2a of the measured central 
velocity dispersion and radial and tangential arc loca- 
tions. Projections of these points in Af c i ustcr -c space and 
M c i ustel -M*/L space are shown in Figure |H The cluster 
Abell 383 does not have any models which lie within la 
of the observations, so the 2a models are plotted instead. 
These projections show that the best-fit models lie on a 
tight relation between M c i uste r and c. From equation ll61 
halos in this mass range should have a concentration be- 
tween 3.2 and 5.0, eliminating most of the best-fit models 
for MS 2137-23 and RXJ-1133. In the case of MS 2137- 
23, this eliminates almost all of the dissipationless mod- 
els (squares). However, weak lensing measurements of 
MS-2137-23 predict a concentration about twice as large 
as simulations, ther eby only eliminating a few models 
(jGavazzi et al.ll200l . 

Also illustrated in Figure 2] is the difference in M*/L 
between the two models. The heavy lines indicate the 
median and 25 — 75 th percentile (SIQR) for the best-fits 
for each toy model. As with the comparison to the L-a 
relation, the dissipational models have lower M^/L val- 
ues than the dissipationless models. For example, the 
median (M*/L) v ratios for MS 2137-23 are 4.3 ±2.4 and 
1.2 ±0.6 for the dissipationless and dissipational models, 
respectively. For RX-J1133, (M*/L) B = 2.1 ±0.8 for the 
dissipational models and (M*IV)b = 5.2 ± 1.7 for the 
dissipationless models. The best-fit dissipational models 
all have an M*/L below 2.5(3.5) in the V(B)-ba.nd. As- 
suming passive evolution, the expect ed value of / L is 
(MJL) B w 4.1 ± 0.95 at z ~ 0.35 dTreu fc Koopmansl 
[2001 Ivan der Wei eTaTl l200i iTreu et al.l I2006D . The 
purely dissipational model for RX-J1133 is therefore 
marginally inconsistent with expected M*/L values. The 
shaded regions in Figure 0] show the expected ranges for 
the stellar mass-to-light ratios for the three clusters; MS 
2137-23 is consistent with both the dissipational and dis- 
sipationless models. Neither the dissipational nor the 
dissipationless models represent an adequate fit to Abell 
383; however, the M*/L values for the dissipational mod- 
els within 2a of the observations are in better agreement 
with the expected M*/L value. 

Taking MS 2137-23 as a specific example, Figure [3] 
plots the stellar and dark matter density of the two mod- 
els. Both models selected have tangential and radial arcs 
and velocity dispersions within the error bars of the ob- 
servations. The best-fit model parameters and the ob- 
served model parameters are given in Table [T] For the 
dissipational model, the dark matter densit y dominates 
over t he stellar density at all radii. From San d et ahl 
(|200l . the best-fit inner slope of the dark matter profile 
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TABLE 1 

Input model parameters for strong lensing 
BCG and cluster: MS 2137-23 



10'° 



Parameter 



Dissipational Dissipationless 



Free (fitted) model parameters 



logM c i uater /M 
concentration (c) 
(M,/L) v 



14.75 
8.99 
1.67 



14.55 
13.27 
5.47 



Observed model parameters 3 



log(L B CG/£©)v- 

^Sersic 
He [kpc] 

(7 [km s" 1 ] 13 



11.5 
4.0 
23.03 
317.33 



11.5 
4.0 
23.03 
321.90 



a Observations of Scrsic fit from [Sand et alj (2004). 
b Observed a = 319 ± 26 km s" 1 dSand et al.1120041) . 
The numbers given are from the best-fit models. 

is 0.57. This is shown along with t he 2a error bars. The 
inner slope derived bv lSand et al.1 (|2004f ) depends on the 
concentration remaining fixed at 400 kpc. If the concen- 
tration is allowed to vary, the best-fit inner slope will 
generally increase by 0.15, bringing it closer to the dis- 
sipationless model. However, the dissipationless model 
also depends strongly on concentration and there exist 
choices for M c i us t cr , M*/i, and c, which fit the obser- 
vations equally well, such that the dissipationless model 
almost matches an NFW profile. 

As illustrated by the strong lensing, the differences be- 
tween the two models for BCGs are small. This is due 
to the fact that the ratio between the stellar mass of 




1 oo 



Fig. 5. — The dark matter and stellar density profiles for two 
best-fit models to the lensing observations of MS 2137-23. The solid 
curve is the dark matter and the lower dotted curve is the stellar 
component of the dissipational model. The dashed curve and the 
upper dotted curve are the dark matter and stars, respectively, of 
the dissipationlcss model. The dash-dotted curve with the shaded 
region is the best-fit model for the dark matter from ISand et alj 
( 2004) along with 2cr error bars. 

the galaxy and the dark matter halo mass is very small, 
on the order of 0.002. Figure [6] shows the ratio of the 
dark matter to stellar mass inside 0.25i? o as a function 
of galaxy luminosity using the scaling relations from H3.1{ 
As the luminosity and cluster mass increase, the differ- 
ences between the two models and the initial NFW profile 
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become smaller. Therefore, the differences between the 
dissipational and dissipationless methods of galaxy for- 
mation will be most pronounced in smaller dark matter 
halos, such as fossil groups and isolated ellipticals. In 
these cases the stellar component of the central galaxy is 
much larger relative to the halo component, making the 
differences between the dissipational and dissipationless 
models more pronounced. 
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Fig. 6. — The ratio of dark to stellar matter contained within 
0.25 Rc as a function of BCG luminosity. The solid line shows 
the dissipational model, the dashed line shows the dissipationless 
model. The dotted line in the center shows the ratio of an NFW 
profile for the dark matter to a Sersic profile for the stars. 



4. COMPARISON TO SAURON DATA 

Figure [6] illustrates that the total mass-to-light ratio 
is an increasing function of galaxy luminosity. This is 
in agree ment with the trend found by the SAURON 
project (jCappellari et al.ll2006f) . which uses integrated- 
field spectroscopic observations of 25 E/S0 galaxies. Us- 
ing these observations and stellar population models to 
determine the stellar mass-to-light rat ios for their sam- 
ple of galaxies, ICappellari et al.l (|2006l ) find that the to- 
tal(dynamical) M/L is consistently larger than M*/L 
and that this difference increases with increasing stellar 
mass. This trend is shown in Figure [7J The ' x '-symbols 
denote the SAURON data and the shaded region is the 
best fit. The stellar mass-to-light ratio (/-band) is never 
larger than 3.4 f or the brightest galaxie s. The best fitting 
line is given by (jCappellari et al.ll2006| ) 



(M/L)j = (2.35 ±0.19) 



0.32±0.06 



(27) 



This fit ignores the galaxy (M32) at Mj ~ -17.5. 

We can compare our dissipational and dissipationless 
models to the SAURON data to determine whether the 
models recover the trend given by equation [27] As in the 
previo us section, we assume t hat each of the 25 galaxies 
in the ICappellari et al.l (|2006l ) study is formed either by 
purely dissipational or purely dissipationless processes. 
We fix the stellar mass-to-light ratio, the total /-band 
luminosity, and the effective radius (R e ) for each galaxy 



Fig. 7. — The trend in M/L with galaxy luminoisty observed 
in the SAURON data (Cappcllari ct al. 20®. The M,/L values 
for these galaxies are all below 3.4. The ' X '-symbols denote the 
SAURON data. The dotted line and surrounding shaded region 
show the best fit to the SAURON data. The triangles (squares) 
show the M/L calculated for the dissipational (dissipationless) 
models for each SAURON galaxy. The solid and dashed lines show 
the best fit for the dissipational and dissipationless models, respec- 
tively. 



to the values given in ICappellari et~aTl ([2006). We then 
vary the dark matter halo mass for the dissipational 
and dissipationless models of each galaxy until the ve- 
locity dispersion within R e for both mode ls matches the 
value reported in ICappellari et all (|2006l ). Since dissi- 
pational formation increases the dark matter content in 
the center of a galaxy compared to the dissipationless 
model, a smaller total halo mass is required to recover 
the same central velocity dispersion in the dissipation- 
ally formed galaxies than in the dissipationlessly formed 
galaxies. The differences in halo mass lead to differences 
in dynamical M/L, which are shown in Figure [7j Both 
the dissipational and dissipationless models reproduce 
the same trend in M/L with luminosity as is shown in 
the SAURON data. However, the dissipationless models 
yield slightly lower M/L values than the dissipational 
models, leading to better agreement with the observed 
values. The standard deviation of the SAURON points 
around the best fit line is 0.11. The standard deviation 
of the dissipationless model points (squares) around the 
best-fit line to the SAURON data (dotted line) is 0.14. 
The same value for the dissipational models is 0.25. How- 
ever, both the dissipational and dissipationless models 
have M/L values higher than those from the SAURON 
data. Therefore, no combination of these models will 
yield the measured SAURON galaxies. However, both 
the dissipational and dissipationless models used here as- 
sume the stellar orbits are isotropic and that the galax- 
ies are spherically symmetric. Both of these assumptions 
will affect the model-calculated M/L values; in the case 
of rotating galaxies, the calculated M/L values will be 
lowered and possibly brought into better agreement with 
the SAURON observations. 

5. EXAMPLE GALAXY: NGC 4494 
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NGC 4494 is an ordinary elliptical galaxy with a B- 
band luminosity of 2.37 x 10 10 L Q . Because it is an iso- 
lated galaxy instead of a BCG, the mass ratio between 
the dark matter halo and the stars is smaller and, there- 
fore, the difference between the dissipationless and dis- 
sipational models of formation will be larger than those 
found for BCGs in massive clusters. As with the BCGs, 
the input model parameters for NGC 4494 can be con- 
strained by observations. 

5.1. Velocity Dispersion using Planetary Nebulae 

Planetary nebulae have been established as a good 
mass tracer in the outer regions of galaxies. These obser- 
vations provide a good compar ison case for our extreme 
models of galaxy formation. iNapolitano et all (2009) 
measure positions and velocities of planetary nebulae out 
to ~ 7 R c in the elliptical NGC 4494, probing the ve- 
locity dispersion for the galaxy much farther out than 
the central velocity dispersion. At these large radii, the 
dark matter will be comparable to, and dominate over, 
the stellar matter (see Figure [9]) , thus emphasizing the 
differences between the dissipational and dissipationless 
models. By fixing the luminosity (Lb = 2.37 x 10 10 L@), 
the effective radius (R c — 3.68 kpc), and the Sersic in- 
dex ( n = 3.30), of the mode l galaxies to observations 
from INapolitano et al.l (|2009f ). we can fit our dissipa- 
tional and dissipationless models to the planetary nebu- 
lae velocity dispersion curves by varying the total halo 
mass and the stellar mass-to-light ratios. These fits 
also include a point for the central velocity dispersion, 
a = 150.2 ± 3.7 km s~\ as r eported in the HyperledcQ 
database (|Paturel et al.ll2003h . The results are shown in 
Figured] For the dissipational model, (M*/L) B = 2.97, 
while the dissipationless model has a mass-to-light ratio 
of 3.96. The total halo masses are 6.0 x 10 11 M for 
the dissipational model and 1.0 x 10 13 M Q for the dissi- 
pationless model. The slightly poorer quality fit for the 
dissipationless model is due to the fact that the galaxy's 
effective radius is close to the minimum allowed for the 
galaxy to form via dissipationless mergers (~ 2.7 kpc), 
as discussed above in §2.41 As the galaxy approaches 
this minimum size, there is insufficient binding energy in 
some of the central dark matter layers to allow the stellar 
layers to cross. The fit of the dissipationless model could 
be improved if we relaxed our model requirement that 
all the energy from the stars is deposited locally, instead 
allowing more energy to be deposited in the outer regions 
of the halo. This could be the case if the orbits of the 
in-falling material were radial orbits instead of perfectly 
circular orbits as assumed in this work. 

Although the dissipational model shown in Figure |8] 
provides a better fit to the data, the -B-band mass-to- 
light ratio required for the dissipational model is signifi- 
cantly lower than the value derived from stella r popula- 
tion models, 4.3 ± 0.7 (|Napolitano et al.ll2009fl . Thus, a 
purely dissipational formation of NGC 4494 appears to 
be ruled out at the ~ 1.9er level. 

Although adjusting the stellar mass-to- light ratio elim- 
inates the differences in the velocity dispersion profile for 
these two models, the difference in dark matter density 
between them remains large (see Figure [9]). The inner 
dark matter density profile for the dissipationless model 

2 http://leda.univ-lyonl.fr 



dissipational, x 2 -1.32, 

a = 150.2, (M./L) B =2.97 _ 
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Fig. 8. — T he velocity dispersion of NGC 4494 from 

INapolit ano et al. (2009) measured using planetary nebulae. The 
two lines show the dissipational and dissipationless galaxy forma- 
tion models. The central velocity dispersion of NGC 4494 is taken 
to be 150.2 ± 3.7 km s" 1 UPaturel et alj|2003ft . 

follows ~ r~°- 2 , while that of the dissipational model 
follows ~ r -1 ' 7 , making the central dark matter densi- 
ties in the two models very different. The slope index 
of the dark matter in the dissipational model is in good 
agreement with that predicted for a final isothermal mass 
distribution in §T] 

5.2. Dark Matter Annihilation 

Although only available in the Milky Way, one 
direct method of probing WIMP dark matter cur- 
rently being explored is the observation of gamma 
rays from the self-annihilation of WIMP dark matter 
particles (iStoehr et al.l 120031: iColafrancesco et al.l [20061: 
iDiemand et al.ll2007D . The signal strength from such an- 
nihilations will be proportional to p^ark- Assuming a 
smooth distribution of dark matter, the ratio of the an- 
nihilation signal strength within an aperture of R c for 
the dissipational and dissipationless models of NGC 4494 
is around 1890, similar to what would be expected for 
a Milky Way sized halo. Although not observable to- 
day, the Fermi gamma-ray space telescope hopes to mea- 
sure the dark matter annihilation signal from our own 
galaxy. The large difference in signal strength between 
the dissipational and the dissipationless models calcu- 
lated here dominates over the boost in signal strength 
expecte d from unresolved s ubstructure, which is of or- 
der 10 (jStrigari et alj|2007t iKuhlen et al.ll2008ft . provid- 
ing another possible test of the formation history of the 
stellar component of galaxies. 

6. CONCLUSIONS 

We have shown that the two extreme cases for the as- 
sembly of the stellar content of galaxies lead to large 
differences in the dark matter density profiles of galax- 
ies, assuming that the initial halo conditions are well- 
described by N-body simulations. The stellar mass den- 
sity dominates over the dark matter density in the cen- 
tral regions of both models; the dark matter density in 
the dissipational models can be as much as two orders 



Dark Matter in Centers of Galaxies 



13 



of magnitude lower at r ~ 1 kpc than the dark matter 
density in the dissipational models. However, because 
galaxies are undoubtedly built up by both dissipational 
and dissipationless accretion, most observations will not 
easily distinguish between these two models. For exam- 
ple, although the best-fit models for BCGs have differ- 
ent stellar mass-to-light ratios, neither is outside the ac- 
ceptable range of values from stellar population mod- 
els and observations. Strong gravitational lensing obser- 
vations of BCGs and their host clusters show that the 
dissipational formation models have M*/L values that 
are marginally too low compared to those expected for 
passively evolving ellipticals. For RX-J1133, the me- 
dian (M*/L) B = 2.1 ± 0.8 and 5.2 ± 1.7 for the dissipa- 
tional and dissipationless models, respectively, while the 
expected (M»/L)b from passive evolution is 4.1 ± 1.0 
(jTreu fe Koopmana 120041 ). This discrepancy between 
Ma,/ L values marginally rules out a purely dissipational 
formation history for BCGs, in agreement with both the- 
oretical expectations and other observational evidence. 
Observations of strong lensing by BCGs have been used 
to effectively rule out dark mat ter density profi les as 
steep and steeper than an NFW (|Sand et al.ll2004l ). fur- 
ther strengthening arguments against a purely dissipa- 
tional formation for the stellar component of BCGs. Al- 
though extreme values for the concentration (~ 10) and 
(M*/L)s (~ 1.0) are allowed by the lensing and dynam- 
ics data, the dissipational model can be ruled out for a 
more constrained and plausible set of model parameters. 

Both models adequately reproduce the trend of in- 
creasing total M/L with galaxy luminosity for E/S0 
galaxies, observed usin g integrated-field spect roscopy by 
the SAURON project (jCappellari et al.ll2006[ ). However, 
the lower M/L values found for the dissipationless mod- 
els are in better agreement with the data. 

Constraints on the stellar mass-to-light ratios can also 
be used to exclude the purely dissipational model of 
galaxy formation in the case of the isolated elliptical, 




r [kpc] 

Fig. 9. — The projected dark matter density profile for the best- 
fit models to NGC 4494. The solid line is for the dissipationless 
model. The long-dashed line is the dissipationless model. The 
short-dashed line is the best fit NFW. The two dotted lines are 
the stellar profiles for the dissipationless (upper) and dissipational 
(lower) models. 



NGC 4494. Fitting the dissipational and dissipation- 
less models to observations of planetary nebulae yields 
(M*/L)s values of 2.97 and 3.96 for the dissipational 
and dissipationless models, respectively. Compared to 
4.3 ± 0.7, the (M*/L)b inferred from stellar synthesis 
models, the purely dissipational model can be ruled out 
at the 1.9<7 level. 

Since the change in the dark matter density for both 
models is directly related to the change in the central 
mass of the halo, the larger the stellar component is rel- 
ative to the dark matter halo, the larger the differences 
between the dissipational and dissipationless extremes 
will be. Therefore, instead of examining the properties 
of BCGs, we propose looking for the differences between 
dissipational and dissipationless formation mechanisms 
using the brightest galaxies of fossil groups and isolated 
elliptical galaxies. The large differences attainable in this 
mass range of galaxies is clearly illustrated by the study 
of NGC 4494. The dark matter density profiles shown 
in Figure [9] have inner slope indices of a « 0.2 and 1.7 
for the dissipationless and dissipational models, respec- 
tively. The differences in the dark matter density pro- 
files for galaxies in this mass range are significant enough 
that they could be probed by galaxy-galaxy weak lensing 
studies, provided the difference in dark matter slopes is 
not removed by averaging over many galaxies with dif- 
ferent formation histories. Finally, the difference in dark 
matter content between the dissipational and dissipation- 
less models yields differences in the signal strength from 
dark matter annihilation of order ~ 1890, far larger than 
the boost factor expected from the unresolved dark mat- 
ter substructure in the Milky Way halo. 

The focus of this paper has been the energetics of the 
dissipational and dissipationless galaxy formation mech- 
anisms, not the mechanisms themselves. For dissipa- 
tional galaxy formation, we have assumed that baryons 
cool and condense in the center of halos, leading to adi- 
abatic contraction of the surrounding dark matter. This 
behavior has been confirmed in cosmological simulations. 
Although simplified, the model presented here is cor- 
rect, on average, for more complicated galaxy forma- 
tion scenarios, including major as well as minor mergers, 
and accretion from filaments instead of spherical shells 
(|Gnedin et al.ll200l . 

The physical mechanism we propose for dissipationless 
galaxy formation is the dynamical friction of small stel- 
lar clumps against a smooth dark matter background. 
In the models used here, we assume circular orbits for 
the incoming stellar material. The inclusion of radial or- 
bits and non-spherical galaxies is left for future work, as 
it requires modeling a combination of dissipational and 
dissipationless formation mechanisms. We assume that 
the build-up of the galaxy occurs via a series of small , 
minor mergers jBezanson et all [2001 iCook et al.1 [2009; 
iNaab et al.ll2009D . not allowing for equal-mass mergers, 
which more violently disrupt the system. Indeed, it has 
been shown in dissipationless N-body simulations that 
equal-mass merger remnants will retain the profile of the 
steepest progenitor dBoylan-Kolchin fc Mdl2004tlDehnenl 
12003 : iKazantzidis et all 120061 : iVass et al.ll2009T ). There- 
fore, cuspy dark matter profiles are robust under major 
mergers. However, if baryons are added to dark mat- 
ter halos, they will presumably condense more than the 
dark matter, making up the bulk of the central, high- 
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density matter in merging halos. As two halos merge, 
the outer, less tightly bound and dark-matter-dominated 
components will be tidally stripped, but the central high 
density, predominantly stellar components will settle into 
the center of the merger remnant, undergoing dynamical 
friction along the way. Therefore, the baryons are an im- 
portant ingredient to dry merger scenarios because they 
ensure the merging clump is sufficiently tightly bound to 
reach the central regions of the nascent galaxy. 

The extreme differences in the inner dark matter halo 
densities for the dissipational and dissipationless mod- 
els emphasize the importance of the addition of baryons 
to dark matter halos. Without introducing modifica- 
tions to the ACDM paradigm, dark matter halo cusps 
can be reduced to cores via the dissipationless forma- 
tion of the central stellar regions of galaxies. The bal- 



ance between dissipational and dissipationless formation 
mechanisms can be probed by observations. Current ob- 
servations of BCGs and ellipticals galaxies are sufficient 
to exclude a purely dissipational formation mechanism 
for these galaxies. Future measurements of stellar mass- 
to-light ratios from microlensing observations, and di- 
rect detection of dark matter in the Milky Way will help 
to constrain the balance between dissipational and dissi- 
pationless formation mechanisms and the dependence of 
this balance on time and environment. 



We thank S. Tremaine and J. Binney for their valuable 
comments and corrections. We would also like to thank 
the referee for his/her comments and suggestions. CNL 
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